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WALL AND BOTTOM HEATING OF LIQUID HYDROGEN IN A PROPELLANT TANK 
by Sidney C. Huntley, James W. Gauntner, and Bernhard H. Anderson 

Lewis Research Center 


SUMMARY 

An experimental investigation was made to determine the behavior of liquid hydrogen 
in a 12 5- gallon scale model of a propellant tank subjected to a range of wall heat flux from 
0. 0012 to 0. 0082 Btu per square foot per second and a range of bottom heat flux from 
0. 0005 to 0. 0098 Btu per square foot per second while discharging from the tank at a rate 
of 0. 04 pound per second under a constant tank pressure of about 2 atmospheres. In- 
creasing temperature stratification in the liquid was encountered with increasing wall heat 
flux; a decrease in stratification was experienced with increasing bottom heat flux. An 
available analysis partially predicted the increase in temperature stratification. 

Generalization of the exit temperature histories showed that the data were applicable 
to other test conditions provided similarity existed in the normalized heat input rate dis- 
tribution. Application of the generalization method to other test data showed that bottom 
heating gave results similar to nuclear heating from below which had a similar overall 
heat input rate distribution. 

No gross changes in liquid behavior were experienced over the range of experimental 
conditions. 


INTRODUCTION 

A knowledge of the effect of heat leaks into the cryogenic propellant in a vehicle 
storage tank is required for the optimization of tank pressurization, tank insulation, and 
residual propellant requirements of the vehicle. In chemical rockets, thermal radiation 
and aerodynamic heating supply a net influx of heat to the propellant from the tank walls. 
In nuclear rockets, additional heating of the fluid arises from the absorption of neutron 
and gamma ray radiation both in the tank walls and in the fluid. This heating of the fluid 
produces both an increase in the bulk temperature and a vertical temperature gradient 
(stratified layer) near the liquid surface. In the case of constant pressure discharge, 
the significance of the stratified layer is realized when a pump is used in the propellant 



feed system. The net positive suction head requirements of the pump limits the permis- 
sible temperature rise in the subcooled liquid. Thus a portion of the stratified layer 
represents unusuable propellant, the quantity of which must be minimized. The amount 
of unusable propellant depends on the heat influx to the tank, both magnitude and distri- 
bution, which in turn is determined by the thermal protection of the vehicle tank and the 
flight environmental conditions. 

Numerous experimental and analytical studies have examined the problems asso- 
ciated with the storage of cryogenic propellants under the conditions of external heat 
loads. Experimental measurements on the transient temperature distribution in the 
liquid are presented in references 1 and 2. These data were obtained under no-flow con- 
ditions. A number of approximate analyses have been developed to predict transient 
temperature profiles in the fluid proper for closed tanks and tanks with outflow. Typical 
examples of this work are presented in references 3 to 5. A visual study of the flow be- 
havior of a noncryogenic fluid in a small two-dimensional vented tank was reported in 
reference 6. The stratification pattern was examined when the liquid heating was from 
the side walls, bottom, or internal absorption of energy or a combination thereof. In- 
frared radiation was used as a heat source. The results indicated that the stratification 
pattern is very dependent on the mode of heating. An excellent review article which dis- 
cusses and summarizes both experimental and analytical work done in this field appears 
in reference 7. A series of tests was conducted under contract to obtain data with wall 
and bottom (including source) heating of liquid hydrogen provided by a nuclear environ- 
ment (ref. 8). Tests were conducted in a 125- gallon tank in which the flow rate, total 
heating rate, and tank pressure were varied. It should be noted that although the total 
heating rate could be varied, there was no provision for separately controlling the side 
wall heating. Presented in reference 9 is a comparison between an analysis and results 
from the nuclear environment tests. The analysis was developed in reference 9 pri- 
marily to predict the transient temperature distribution in liquid hydrogen with wall and 
nuclear heating for conditions of constant pressure discharge. 

To further study the liquid behavior over a wider range of heating rate distribution, 
a series of experiments was conducted at the Lewis Research Center with liquid hydrogen 
in a 125-gallon tank similar to that used in the nuclear environment tests. This series 
of tests used separately controlled radiant heaters located along the side and bottom of 
the tank to permit variation of both side and bottom heating. It is the purpose of this re- 
port to document the results of this series of tests. The analysis presented in refer- 
ence 9, giving general agreement with the experimental nuclear data, was applied to the 
present series of tests and its range of applicability is briefly discussed. A set of non- 
dimensional scaling parameters are developed and utilized in the presentation of the ex- 
perimental data so that they can be made more directly applicable to other test condi- 
tions. The generalization was then applied to the results of the nuclear environmental 
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tests and comparisons are made with the data reported herein. 


APPARATUS AND PROCEDURE 

The experiment was designed to simulate a rocket vehicle system. The experi- 
mental tank and test parameters were chosen by scaling a rocket vehicle tank having a 
capacity of 100 000 pounds of liquid hydrogen. The value of the flow rate was chosen to 
yield about the same time to outflow as the full size tank and the tank was pressurized to 
about 2 atmospheres, which was deemed adequate for typical missions. The test tank 
used in this series of tests was geometrically the same as the one used in the nuclear 
environmental tests (ref. 8). This was deemed necessary so that a direct comparison 
could be made between the results reported herein and those reported in reference 8. 

The magnitude of the heating rate applied to the tank was chosen to preserve heating 
loads (heating rate per unit volume of liquid) for typical chemical and nuclear rocket 
missions for the full-size tank. (Symbols used in this report are listed in appendix A. ) 

The experimental apparatus is reported in detail in appendix B; thus, only a cursory 
description is included here. The test apparatus (fig. 1) consisted of a test tank enclosed 
in a vacuum chamber with provisions for supplying liquid hydrogen, pressurizing and 
venting the tank, throttling the outflow, and providing separately controlled heating of 
the tank side wall and bottom. A vacuum system was provided that had the capability of 

c 

providing a vacuum of less than 10”^ torr, thereby reducing the conductive heat leak 
through the vacuum space to a negligible amount. The outer shell of the vacuum space 
was provided with liquid nitrogen cooling to reduce the radiant heat flow into the sides 
and bottom of the test tank. Solar radiation and/or aerodynamic heating was simulated 
by using electrical radiant heaters (fig. 2) suspended in the vacuum. A detailed descrip- 
tion of the radiant heaters is also presented in appendix B. 

Two temperature measurement systems were used for the tests: one system con- 
sisted of carbon resistor thermometers to measure temperatures in the cryogenic range 
and the other of copper constantan thermocouples to measure temperatures above the 
cryogenic level. Temperature data are presented in this report for thermometers 
located at 0, 16. 9, 24. 7, 32. 7, and 40. 7 inches above the exit port along the tank center 
axis. Other types of measurements included tank pressure, liquid flow rate, and LLiiuid 
level position. 

The flow test conditions for each heating rate distribution considered are presented 
in table I. Each test is characterized by a different heating rate distribution which is 
identified by the average wall heat flux relative to the average bottom heat flux; that is, 
the identification H - L indicates that the average wall heat flux (H) was somewhat 
higher than the average bottom heat flxax (L). Intermediate values of heat flux are 
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designated as either ML or MH. The initial height of liquid, initial liquid temperature, 
tank pressure with resultii^ saturation temperature, and flow rate are also shown. 

Three types of tests were performed in this investigation: boiloff, no-flow, and flow 
tests. Detailed operating procedures for each of these tests are presented in appendix B. 
The boiloff and no-flow (self-pressurizing) tests were used as calibration runs to sub- 
stantiate the heating rate distribution for the flow tests. 

The method of obtaining the heat input distribution for the flow tests was based on a 
calculation of the net radiant heat exchange between the radiant heaters and the tank 
walls. The net heat input rate from the heaters was obtained from the equation 





q(x) ^ dx 


( 1 ) 


The quantity q(x) is the net heat flux from the radiant heaters (see eq. (Cl) of appen- 
dix C), S is the wetted surface area of the tank, and x^ is the location of the liquid sur- 
face. Because the calculations did not consider other forms of heat transfer that were 
present (e. g. , radiation from the tank dome), independent checks of the heat input were 
performed and corrections to the heating rate were made. This was accomplished by 
first integrating the radiant heating rate distribution over the time to outflow for each 
run. The average heat input rate was established by means of the ejq>ression 


^H“ t 

m 



dS 

q(x) ~ dx dt 
dx 


( 2 ) 


The average heating rate of the liquid for each run was obtained by the integration of 
the temperature history at the exit port using the equation 


pc V 

\=~f~ i ^T(o,t)dt (3) 

m 0 


where the density p and specific heat at constant pressure Cp were evaluated at the 
average temperature over the time to outflow defined by the ejq)ression 


T = 


T. + “ 
in ^ t 


m 



AT(o, t) dt 
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and ^ is the average volumetric flow rate calculated from the slope of a straight-line 
cxirve fit of liquid volume history. The liquid volume history was obtained by observing 
the time at which the liquid surface passed the several temperature sensor locations. 

The average heating rate from the heaters Qjj and the average heating rate of the 
the liquid thus obtained were then used as a basis for adjusting the calculated radiant 
heat distribution by means of the equation 


♦ 

% 


( 4 ) 


» “ 

The value of /^jj obtained for each flow test is shown in table I. In order to estab- 
lish the validity of this adjustment, independent checks of the total heating rate 
and the heating rate distribution were made using boiloff and no-flow tests. A 

detailed discussion of the heating calculations and the methods of substantiating these 
calculations appears in appendix C. 

The heating rate distribution (eq. (4)) and the flow test conditions (table I) for each 
heating configuration were used as input for the analysis developed in reference 9 and 
the liquid temperature history was calculated. A detailed discussion of the method of 
application and the basic equations used in the analysis is presented in appendix D. 


TEST RESULTS 

Several flow tests were conducted in which the side wall and bottom heat flux were 
held constant. Each test had a different heating configuration distinguished by the rel- 
ative values of wall and bottom heat flux (table I). The tests were conducted primarily 
to show the effect of varying wall or bottom heat flux on the behavior of liquid hydrogen 
during constant pressure discharge. Heat input rate distributions obtained from the 
flow test experiments are given, followed by the presentation of the experimental tem- 
perature histories and a comparison between the experimental results and anal 3 rtical 
predictions. 


Heat Input Distribution 

The heat input rate from the heaters is presented in figure 3 as a function of liquid 
height for the several configurations. Data points represent calculated values of heat 
input rate using equation (1). Solid curves are used to depict test configurations related 
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by having increasing values of wall heat flux, while the dashed curves depict increasing 
values of bottom heat flux. The lowest curve represents the condition of minimum heat 
leak to the tank obtained with no power supplied to the heaters; the top curve was ob - 
tained with both heaters giving about equal heat fluxes, such that the bulk of liquid was 
heated to saturation by the time the tank was emptied. 

Time-averaged heating rates in the liquid and from the heaters are presented in 
figure 4 for the several heating configurations. The quantity - Qjj)/Qjj defined as 
the heating error is shown as a function of the average heating rate from the heaters. 

Also shown are the results obtained from the no-flow tests described in appendix C. The 

-* • 

dashed curves represent constant - Qjj values of 15 and 30 watts. The data for the 
two types of tests show good general agreement although the heating error was greater 
for the no-flow type of tests. 


Experimental Temperature Histories 

The temperature history at several axial locations in the tank is presented in fig- 
ure 5 for each heating configuration (see table I for initial conditions). The experi- 
mental data are depicted by solid curves; the dashed curves, which show predicted 
temperature histories, will be discussed in the next section. In general, the tempera- 
ture at a given location rises gradually with time during the outflow of liquid from the 
tank. The rapid increase in temperature indicates the presence of a stratified layer of 
warm liquid near the surface. The nearly constant temperature after the rapid increase 
corresponds to saturation temperature. Apparently the dip in the temperature history 
occurs when the thermometer becomes exposed to the gas and evaporization of the liquid 
film causes a slight cooling effect. 

The data for the L - L configuration which represents the minimum heat input is 
presented in figure 5(a). Under this heat input condition, the gain in liquid temperature 
at the exit was small until the tank was nearly empty. The constant temperature at any 
time for several locations indicates that the heat was distributed uniformly throughout 
the liquid except near the surface. 

The H - H configuration (fig. 5(b)) had the maximum heat input to the tank and was 

chosen to result in heating the bulk of liquid to saturation temperature just prior to the 

tank being empty. Comparison of the results of this test with those of the minimum heat 

-3 

configuration (fig. 5(a)) shows that the increase in average heat flux from 2. 4x10 
to 8X10' Btu/per square foot per second resulted in a general increase in liquid tem- 
perature at any given time. A measurable axial temperature gradient existed in the 
liquid with the H - H configuration as is shown by the spread in temperature of several 
thermometers at a given time. 
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Comparison of the data for the L - L, ML - L, MH - L, and H - L configurations 
(figs. 5(a), (c), (d), and (e)) shows the effect of increasing the wall heat flux from ap- 
proximately 2. 8x10 to 7x10“ Btu per square foot per second with low bottom heat flux. 
A comparison of the data shows that the exit temperature rise in about the first 300 sec- 
onds was approximately the same although wall heat flux was increasing. Since the in- 
crease in wall heat flux occurred above a height of about 16 inches from the tank exit 
(fig. 3), the first 3 cubic feet of liquid to leave the exit would not sense the increase in 
heat. Therefore, with an outflow rate of about 0. 04 poimd per second, about 300 sec- 
onds should be required before the influence of wall heat is observed at the exit. In- 
creasing the wall heat flux resulted in an increasing rate of exit temperature rise (after 
300 sec) except when the tank was nearly empty. The spread in temperature between 
given locations in the tank increased with increasing wall heat flux at a given time indi- 
cating more heat is accumulating in the stratified layer. 

A comparison of the data for the L - ML, L - MH, and L - H configurations 

(figs. 5(f), (g), and (h)) shows the effect of increasing the bottom heat flux from about 
-3 - 3 

3. 6x10 to 9. 8x10 Btu per square foot per second with low wall heat flux. As in the 
case of variable wall heat flux, increasing the bottom heat flux results in an increasing 
rate of exit temperature rise but, since the heat is added at the tank bottom, the spread 
in temperature between thermometers is smaller. The liquid temperature at each loca- 
tion rises at a nearly uniform rate with bottom heating. Apparently incr easily the heat 
at the tank bottom results in increased mixing with the fluid approaching a completely 
mixed condition except for a thin layer near the liquid surface. 


Comparison of Experimental and Predicted Temperature 

Predicted temperature histories for each of the experimental heating configurations 
based on the analysis of reference 9 (see appendix D) are in general agreement with the 
experimental data (fig. 5). Deviations that occur may be classified into two types: one 
pertaining to mass and the other to energy. Deviations due to mass occur primarily be- 
cause of the difference between using an average flow rate in the analysis and the in- 
stantaneous flow rates of the experiment. Deviation of the actual from the predicted 
time for the surface to pass any location in the tank is consequently due to flow rate 
variations from the average value during a test or to slight differences between actual 
probe locations and those used in the analysis. 

Deviations pertaining to energy refer to the suitability of the analysis to predict the 
liquid behavior. The predicted temperature history at the tank exit, based on considera- 
tion of energy conservation, should prescribe the total input energy to the system. The 
determination of the input energy, on the other hand, is dependent on both the time- 
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averaged experimental exit temperature history and the time- averaged volume history. 
Temperature deviations at the tank exit therefore result from both flow rate variations 
and the manner in which the analysis distributes the energy input at any time. 

Slight deviations pertaining to mass between the experimental data and the analysis 
result in large temperature deviations because of the rapid temperature rise that occurs 
as the liquid surface approaches a given sensor location. Temperature deviations of 
this type appear to be of some consequence but are not considered to be of major signi- 
ficance because refinements of the analytical technique are possible. Of greater interest 
here are the temperature deviations, between the experimental data and the analysis, 
pertaining to energy. The analysis prescribes a temperature profile dependent on heat 
flux distribution along the wetted tank surface and a percentage of the total heating rate 
which enters the stratified layer. The temperature profile thus established determines 
the predicted temperature rise at any given location. 

A comparison of the predicted and experimental temperature history at several loca- 
tions for the cases with increasing wall heat flux (L - L, ML - L, MH - L, and H - L, 
see figs. 5(a), (c), (d), and (e)) shows that deviations in temperature profiles increase 
with increasing wall heat flux. A similar comparison between cases with increasing 
bottom heat flux (figs. 5(f), (g), and (h)) indicates the opposite trend. It is apparent that 
best agreement between the analysis and data is obtained with heat input distributions 
having relatively high bottom heating. However, the general agreement between the pre- 
dicted and measured temperature histories indicate the analysis of reference 9 may be 
useful for interpolation. 


GENERALIZATION OF DATA 

For experimental data obtained in small tanks to be applicable to a full-scale tank, 
such tests must simulate the operating conditions prevailing in the environment in which 
the full-scale system will function. This requirement necessitates the establishment of 
a set of dimensionless parameters that must have the same value for both the model and 
the full-scale system. Such a set of dimensionless parameters is derived to obtain 
generalized temperature histories. 


Dimensionless Parameters 

Dimensionless parameters are desired to describe the system consisting of the 
liquid bounded by the tank walls and the liquid gas interface at any time during flow at 
constant pressure. An energy balance of the system may be expressed in terms of the 
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time rate of increase in liquid enthalpy, the rate at which enthalpy is being transported 
from the system, and the heating rate entering the system. When no heat transfer at the 
mterface and no radial temperature gradients are assumed, the energy equation can then 
be written as 




where w is the flow rate and V(x_) is the volume of liquid in the tank at any given time. 

• s 

A saturation heating rate Q_ is now defined as the energy required to heat the initial 

s 

liquid content to saturation temperature in the same time required to empty the tank. 
This saturation heating rate can be illustrated by means of the expression 


pcV(L)AT 

'*8 = ■ V 

m 


( 6 ) 


Combining equations (5) and (6) yields the energy equation as 



dv + (o, t) = ^ — 


. Qfi (x„) 


Q«(L) 


( 7 ) 


where 


T 



V(L) 
i> = ^ 
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In this form the energy equation states that the nondimensional temperature history 
everywhere in the tank is a function of the nondimensional heat input. rate distribution 
(Xs)/^jj (L) and a heating parameter The heating parameter ^ relates the initial 
heating rate to the other test conditions of flow rate w and tank pressure (in terms of 
saturation temperature AT^). 


Heat Input Rate Distribution 

The heat input rates to the liquid, normalized to the initial values of heat input rate, 
are presented in figure 6 for the several heating configurations. Also shown by the 
dashed line is a typical total heating distribution curve from the nuclear heating data of 
reference 8 (this curve will be discussed in a later section). 

The generalized heat input rate distribution for the H - H configuration, having been 
obtained with about constant heat flux, represents the nondimensional tank surface area 
as a function of x/L (see eq. (1)). Consequently, curves which lie above represent con- 
figurations with relatively higher bottom heat flux; those below represent relatively 
higher wall heat flux. 

Increasing the bottom heat flux (L - ML, L - MH, L - H) from 3. 6x10"^ to 9. 8x10“^ 
Btu per square foot per second did not result in as large a change in shape of the gen- 
eralized heat input rate distribution curve as expected (fig. 6). This was partially due 
to an increasing wall heat flux with increased bottom heat flux. The heat input rate dis- 
tribution curves have about the same general shape as the total heating distribution of the 
nuclear data, which is shown on figure 6 but will be discussed in a later section. 

An increase in wall heat flux (L - L, ML - L, MH - L, H - L) from 2. 8x10"^ 

O 

to 7x10" Btu per square foot per second also did not affect the general shape of the dis- 
tribution curves. In general, the heat input rate distribution curves can be classified 
into two general shapes, those with high bottom heating and those with high wall heating. 
This general groupii^ of the data is used as a convenient basis for the following discus- 
sion. 


Generalized Exit Temperature Histories 

The input to the generalized temperature histories was considered to consist of two 

•' • « 

parts, 2 and Qj^(Xg)/Qj^(L). The latter, discussed in the previous section, was shown 
to consist primarily of two types of heating, either high bottom or high wall heating. The 
generalized exit temperature histories for each type of heating are presented in figure 7. 
The heating parameter 2 is shown for each case. Increasing the initial heating rate to 
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the system with the same flow rate and tank pressure results in higher values of the 
heating parameter ^ . As expected, increasing this parameter increased the rate of 
exit temperature rise for either type of heat input distribution. With high bottom heating 
a greater portion of the initial heating rate enters the liquid for a longer period of time 
than with high wall heating (fig. 6). Consequently, exit temperatures with high bottom 
heating (fig. 7(a)) are greater than those with high wall heating (fig. 7(b)) at any given 
time for about the same heating parameter 

Neither the wall nor the bottom heating data (fig. 7) are strictly families of curves 
since the heating rate distributions differ slightly with each type of heating (fig. 6) in 
addition to variations in the heating parameter. For this reason, caution should be ex- 
ercised in interpolating these data to intermediate values of heating parameter. 


Comparison of Bottom and Nuclear Heating Data 

A typical generalized heating distribution curve from the nuclear heating data of ref- 
erence 8 was presented in figure 6. The shape of this curve is similar to those of the 
present data having high bottom heating. It was considered of interest, therefore, to 
compare the present bottom heating data with the nuclear heating data. The nuclear heat- 
ing rate distribution was essentially independent of initial heating rate, and consequently 
the generalized total heating rate distributions all had about the same shape. Therefore, 
the exit temperature histories should form a consistent family of curves with the heating 
parameter ^ , 

The generalized exit temperature histories of the nuclear heating data are presented 
in figure 8. Also shown are the corresponding operating conditions for each test rvin and 
the heating parameter The several test runs have been presented in the order of de- 
creasing values of the heating parameter. Four test runs yielded heating parameters 
from 0. 373 to 0. 388, which covered a range of two in pressure level and in initial total 
heating rate. In another instance, heating parameters of 0. 117 and 0. 120 occurred with 
the effect of initial total heating rate variations approximately balanced by changes in tank 
pressure and flow rate. In both instances, the generalized exit temperature histories 
were consistent showing that the heating parameter adequately accounted for changes in 
operating conditions. 

A comparison of the bottom heating data with the nuclear data (figs. 7(a), and (8)) 
shows consistent generalized exit temperature histories with the parameter Q, . The 
L - ML configuration of figure 7(a) had a heat input rate distribution most similar to the 
nuclear data. In this instance, the shape of the generalized exit temperature history 
(^ = 0. 117) bears a striking resemblance to a nuclear data run having the same value of 

The same value of 2 was achieved in the nuclear data run with a greater initial 
heat input rate but with a compensating higher flow rate and tank pressure. Perhaps the 
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most significant factor of the comparison of bottom and nuclear heating is that the ap- 
proach of simulating nuclear heating by bottom heating appears to be feasible. 


SUMMARY OF RESULTS 

The following results were obtained from an experimental study of the behavior of 
liquid hydrogen contained in a tank subjected to wall and bottom heating during flow at 
constant pressure: 

1. Increasing the wall heat flux with low bottom heat flux showed that the added heat 
accumulated in a stratified layer. A decrease in stratification occurred when the bottom 
heat flux was increased with a trend toward a completely mixed condition except for a 
thin surface layer. 

2. No gross changes in liquid behavior were experienced over the range of experi- 
mental conditions. 

3. An available analysis partially predicted the increase in temperature stratification 
with increasing wall heat flux. The general agreement between the predicted and 
measured temperatures may be useful to interpolate between experimental results. 

4. A set of dimensionless parameters was derived to obtain generalized exit temper- 
ature histories. Use of the parameters showed that the data can be made applicable to 
other test conditions provided similar heat input rate configurations are used. 

5. Application of the generalization method to other test data showed that bottom 
heating gave exit temperature histories similar to nuclear heating from below with 
similar overall heat input rate distribution. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 23, 1965. 
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APPENDIX A 


A tank cross-sectional area 

c specific heat at constant pressure 

E interchange factor, liquid energy 

e energy per unit mass 

f(t) ratio of exit temperature rise to 
saturation temperature rise 

hy heat of vaporization 

L initial height of liquid 

n exponent, eq. (D2) 

heating rate per unit volume 

tank pressure 

Q heating rate 

^ heating parameter, ratio of initial 

heating rate to saturation heating 
rate, eq. (7) 

Qjj heating rate from heaters, eq. (1) 

heating rate of liquid, eq. (3) 

Q_ saturation heating rate, eq. (6) 
s 

q heat flux 

r radius 

r^ radius of tank bottom 

S surface area of tank 

Sj surface area of tank from 

X = 13. 5 to 

T temperature 

t time 

V volume 


volumetric flow rate 

ratio of volume to initial volume 
of fluid 

w weight flow rate 

X axial distance from tank bottom 

6 thickness of stratified layer 

e emissivity 

0 ratio of temperature difference to 

saturation temperature dif- 
ference 

p fluid density 

0 Stefan- Boltzmann constant 

T ratio of time to time to outflow, 

!p similarity parameter, eq. (D2) 

Subscripts: 

1, 2, 3 tangent points 
H heater 

1 bottom of stratified layer 

in initial condition 

I liquid 

m maximum 

s saturation, surface condition 

T tank 

w wall 

si stratified layer 

Superscript: 

~ average conditions 


SYMBOLS 

V 
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APPENDIX B 


DESCRIPTION OF TEST APPARATUS, INSTRUMENTATION, AND PROCEDURE 

Tank 

The tank, constructed of 304 stainless steel, had a liquid capacity of about 125 gal- 
lons. The tank geometry consisted of a 32-inch-diameter cylinder mounted on top of a 
spherical zone transition to the frustrum of a 45° half angle cone on a spherical segment 
(see table n). A 3/4-inch schedvile 40 liquid outflow pipe extended below the tank. The 
tank skin was 1/8 inch in thickness. A flange and transition section were located at the 
upper end of the cylinder to permit removal of the tank from the hemispherical top. The 
outer surface of the tank (and the inner simface of both the wall and bottom heaters) was 
sand blasted and spray painted with a flat black lacquer to provide a high emissivity sur- 
face. 

A window-type opening was available in the side of the tank, but it was sealed with a 
steel cover for this set of experiments (fig. 2). Two similar ports installed in the tank 
dome were used for the purposes of viewing and illumination of the tank interior. 


Radiant Heaters 

The radiant heaters were constructed from 304 stainless steel ribbons. The wall 
heater consisted of a double lead coil of 1- by 1/16-inch ribbon cylindrically wrapped 

with a 37^ - inch inside diameter and a 2^ - inch pitch. Each coil had approximately 

13 convolutions. The bottom ends of each coil were joined with a 2^- inch strip of the 

heater material to form one continuous strip to which electrical power leads were at- 
tached at each end. A 13-inch-diameter opening had to be provided for the protruding 
flange on the side of the tank. Sections of the ribbons were removed to provide the open- 
ing, but the concept of a continuous strip was maintained by rejoinii^ appropriate bands 
(see fig. 2). The overall height of the wall heater was about 30 inches. Integrity of the 
coil was maintained by using vertical insulating bars riveted to each ribbon. The bars 
also provided a means of supporting the heaters from the top of the vacuum jacket while 
maintaining both electrical and thermal barriers. 

The bottom heater consisted of a double lead coil of 3/4 by 1/16 inch ribbon helically 

wrapped with a 45° half angle from the axis with a major diameter of 36^ inches using a 
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1— - inch pitch. Each coil had approximately 13 convolutions with the bottom lead ends 

4 

joined to form one continuous strip. The bottom heater was suspended from the wall 
heater by using similar bars for support and integrity as in the case of the wall heater. 

The double lead coils of one continuous strip were used for both the wall and bottom 
heaters to prevent inductive coupling between the heaters and the tank. Voltage regu- 
lated alternating current was separately supplied to each heater. 


Instrumentation 

Two temperature measurement systems were used for the test; one consisted of 
carbon resistors to measure temperatures in the cryogenic temperature range and the 
other of copper constantan thermocouples to measure temperature above cryogenic tem- 
peratures. The temperature measurement positions are shown in table III. Thirty 
carbon- resistor thermometers were mounted on a rake extending down the tank axis. 
These thermometers were spaced at 1/4-inch increments near the liquid surface of the 
full tank at the 4-inch increments otherwise. A thermometer placed at the tank exit was 
used to observe the liquid temperature leaving the tank. Eight additional thermometers 
were mounted radially 1/8 inch and 1/4 inch off the tank wall. The 100-ohm, 1/10-watt 
resistors were mounted in a horizontal position. Thirty-six copper constantan thermo- 
couples, mounted along the entire axis of the tank at 4 inch or less increments on the 
axial rake, were used to observe gas temperatures. The thermocouples were spaced at 
l/4-inch increments near the surface of the full tank and near the surface of the half-full 
tank. Four additional thermocouples were mounted radially 1/8 inch and 1/4 inch off the 
tank wall. Six thermocouples were attached to the outer wall of the test tank to measure 
tank wall temperatures. Three thermocouples were attached to the inner surface of each 
radiant heater to measure the heater temperature. 

The response time of both the carbon resistor thermometers and the thermocouples 
was expected to be at least 2 seconds. Previous experience with this type of carbon re- 
sistor gave indication that an accuracy within 0. 1° R was obtainable in the range of liquid 
hydrogen temperatures. The thermocouple accuracy was estimated to be within 1° R at 
temperature levels greater than 138° R. (An atmospheric liquid nitrogen reference bath 
was used for the thermocouples. ) 

A strain-gage type pressure transducer was used to measure tank pressure. A 
venturi, installed in a vacuum chamber several pipe diameters downstream of the tank 
shutoff valve, together with an upstream pressure transducer, a differential pressure 
transducer, and a downstream carbon resistor thermometer were used to indicate the 
liquid flow rate. 
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Test Procedure 


The general procedure of the investigation was to first establish a heating configura- 
tion consisting of a particular wall heat flux and a particular bottom heat flux into the 
test tank. With each configuration, a flow test was conducted which consisted of filling 
the tank with liquid hydrogen, pressurizing to about 2 atmospheres, and then maintain- 
ing the pressure constant during outflow at a rate of about 0. 04 pound per second. The 
flow test conditions for each configuration are presented in table I. 

Boiloff and no-flow tests were performed to check the method of obtaining the heat 
input rate distribution. Boiloff tests were performed in which the tank, filled with liquid 
hydrogen, was allowed to vent at a constant tank pressure. During these tests the liquid 
level history was obtained. 

No-flow tests were performed which consisted of establishing a given heating con- 
figuration with the tank vented after which the fill and vent valves (see fig. 1) were 
closed; this allowed the tank to pressurize. The liquid temperature and tank pressure 
history were used to determine the heat input rate associated with a particular mass of 
liquid. 

The double wall of the outer shell (see fig. 1) was filled with liquid nitrogen to en- 
sure a constant reproducible wall temperature which would yield a low ambient heat 
leak. To decrease the heat leak from the flanges, a tub-like structure was built around 
the exit port and filled with liquid nitrogen. 

Prior to the start of the flow tests, the throttle valve was opened and liquid hydrogen 
from the supply dewar was flowed through to chill the venturi for approximately 15 min- 
utes prior to start of outflow from the tank. At this time, the throttle valve was ad- 
justed to yield the desired flow rate. About 20 seconds prior to start of flow, the tank 
was pressurized with hydrogen gas to 1 atmosphere above the initial pressure level. 
(Prior to tank pressurization, the liquid contents were at saturation temperature cor- 
responding to the initial pressure. ) Next, the fill valve was closed. A few seconds after 
the tank had achieved stable pressure, the shutoff valve was opened and the flow of 
hydrogen began. Additional pressurizing gas was used as needed to maintain a constant 
tank pressure during the entire flow run. The pressurant was introduced into the tank 
through a gas diffuser to avoid disturbing the liquid surface. 
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APPENDIX C 


HEAT INPUT RATE DISTRIBUTION 
Heat Input From Radiant Heaters 

Preliminary runs with the test apparatus indicated a problem not previously antici- 
pated. Cooldown of the radiant heaters resulted in a temperature gradient in each heater 
in the vertical direction. The low value of thermal conductivity of the 304 stainless steel 
heater, combined with the long length of each coil, required an excessive time period to 
equalize the temperature throughout either heater. Therefore, a simple thermal radia- 
tion calculation was used to account for the variation in heating rate resulting from the 
temperature differences. 

Each heater was considered to consist of three zones, each zone having a uniform 
temperature measured by the thermocouple in the zone. The net radiant heat exchange 
was then calculated between each zone of a heater and the corresponding tank surface, 
assuming that both the heater and the wall were gray surfaces and the close spacing of 
the heater to the tank wall was analogous to concentric cylinders. The net heat exchange 
per unit tank area per unit time from each zone was then computed from the equation 

q = ctE(T^ - T^) (Cl) 

where a is the Stefan- Boltzmann constant and 



Constant emissivities of 0. 8 were assumed for both surfaces. The net heat input 
rate to the liquid from the heaters was then obtained from the equation 





(C2) 


where q is obtained from equation (Cl) and x is the height of liquid in the tank. 

o 
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Boiloff Tests 


A boiloff test was conducted in which the liquid level history was measured during 
the time required to completely boil off the initial tank contents. The heat input rate 
from the heaters was also calculated for this test in the manner previously described 
(eq. (C2)) by using the measured heater temperatures. This heat input rate from the 
heaters was used to calculate the time required to boil off a given height of liquid. The 
equation used for this calculation was 


t = 



h^A(Xg) 

Qh<==s> 


dxg 


(C3) 


The heat input rate to the liquid during the boiloff test was calculated from the measured 
liquid level history using the rate of change of liquid height in the expression 


dx 

% = dT 


(C4) 


The experimental results of the boiloff test are presented in figure 9 as the history 
of liquid height. Also shown is the calculated history of liquid height based on the heat 
input from the radiant heaters obtained from equation (C3). A comparison between the 
two curves indicates the method used to calculate the heat input rate from the radiant 
heaters, when applied in this manner, resulted in predicting within 5 percent the time re- 
quired to completely boil off the initial tank contents. The heat input rate distribution 
calculated from the radiant heater method (eq. (C2)) is presented in figure 10 as a func- 
tion of liquid height. Also shown is the calculated heat input rate to the liquid obtained 
from using the observed history of liquid height (eq. (C4)). The two methods of calculat- 
ing heat input distribution were considered to agree satisfactorily within the accuracy of 
measurement. At a height of 45 inches, the liquid heating rate showed about 10-percent 
higher power level than indicated by the heater calculation. This test was conducted at a 
high power level (1 kW into the full tank) to keep the test duration within a reasonable 
time period. 


No- Flow Tests 

No-flow tests were conducted to check the method of calculating the heat input rate 
from the radiant heaters. The test conditions for the no-flow tests are summarized in 
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table IV. In a few instances, two tests were conducted using the same heating configura- 
tion but with different initial heights of liquid to check heat distribution. 

The energy of the liquid at a given time t during a no-flow (or self-pressurizing) 
test was calculated from the equation 

E(t) = /* ® p(x,t)e(x,t)A(x)dx 
*'0 


where the local density p(x, t) and specific energy e(x, t) were obtained from National 
Bureau of Standards literature as a function of tank pressure and a local axial liquid 
temperature T(x, t). At a given time, the temperature data was curve-fit using a three- 
point, one-slope fit for each region defined by thermometer locations. An average heat- 
ing rate of the liquid (third column of table IV) was then calculated from the equation 



E(t J - E(o) 


The heater temperatures deviated a small amount during the course of a test. An 
average heat input rate from the radiant heaters (fourth column of table IV) was there- 
fore used to account for this slight deviation. The equation used to obtain this average 
heat input rate was 


. 1 /•'■» A'*' dS 

where q(x) was obtained from equation (Cl) and the time period of the no-flow test t^ 
(fifth column of table IV). 

The average heat input rate from Ihe heaters was consistently less than the average 
heating rate of the liquid (table IV). The heating error defined by the ratio (Q^ - Qjj)/Qh 
is presented in figure 11 as a function of the average heat input rate Qjj. Also shown is 
a dashed curve representing a constant difference of 30 watts. The dashed curve repre- 
sents a mean average of the experimental data within the accuracy of measurement. The 
general consistency between the heating errors for both heights of liquid with the same 
configuration was considered sufficient evidence to warrant using the radiant heat calcu- 
lation method as a means of determining the heat input rate distribution. 
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APPENDIX D 


APPLICATION OF ANALYSIS 

A method for the anal 3 dical prediction of liquid temperature history subjected to 
heating during flow from the propellant tank at a constant flow rate and pressure was 
developed in reference 9. The use of the method requires a knowledge of tank geometry, 
test conditions, liquid properties, heat input distribution, and assumptions regarding the 
manner in which the heat input will influence the liquid behavior. Some of the basic equa- 
tions used in reference 9 are presented herein for completeness, and the application of 
the analysis to the present data is presented. The system was considered to include only 
the liquid in the tank while discharging at constant pressure and flow rate. The boundary 
conditions included heat transfer from the walls, the surface temperature (corresponding 
to tank pressure), and heat added to the system by the absorption of nuclear radiation. 

The anal 3 dical flow model consisted of a lower region of completely mixed liquid with an 
upper region of stratified liquid having a temperature gradient extending from the com- 
pletely mixed liquid temperature to saturation temperature at the surface. 

The basis for the flow model of the analysis was to treat the nuclear heating in the 
liquid and a portion of the heat transfer from the bottom section of the tank as though it 
resulted in completely mixing the liquid in the region below the stratified layer. The 
heat transfer along the vertical walls was then considered to form the stratified layer, 
which was assumed to have the property of similarity (that is, the property that two 
temperature profiles AT(x, t) at different times differ only in a scale factor in x and 
AT). By this assumption, the temperature in the stratified layer can be written in the 
form 


AT(x, t) = ATg |f(t) [l - ^(x, t)j + tpix, t)| (Dl) 

where ^(x, t) is the similarity parameter and the parameter f(t) j^l - tp{x, t)J can be 
interpreted as the necessary contribution of nuclear (and bottom) heating in the strati- 
fied layer to preserve the condition of similarity. By specifying that the parameter ip 
identically vanishes below the stratified layer, equation (Dl) can be used as the dimen- 
sionless profile for the entire fluid. It is apparent that ip must also satisfy the condi- 
tion that at the liquid surface ^ = 1. It was foimd in reference 9 that under the condi- 
tions where the side wall heating was small relative to nuclear (and bottom) heating, the 
parameter ip(x, t) may be approximated by the form 


ip{x, t) = 



(D2) 
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where x.(t) is the lower extremity of the stratified layer and 6(t) is the thickness of the 
stratified layer. To determine the three unknowns, f(t), 6(t), and the exponent n, the 
energy equation of the system and the assumed energy balance between the boundary and 
the wall heating of the forms 

were used. The assumption was made that the contributions of side wall heating and 
nuclear heating (including bottom heating) could be micoupled. When equation (Dl) is 
substituted into (D3), the uncoupling results in 

pCp ATg ^ i(t)(l - !f)dV + Cpi ATgf(t)=y" PpdV (D5) 

^“P 1 4 * ' 4 

The system of equations (D4), (D5), and (D6) was then solved for the unknown quantities 
f(t), 6(t), and the exponent n. 

For the comparison of analysis and experimental data presented herein, the volume 
integral on the right hand side of equation (D5) was replaced by a surface integral of heat 
flux over the tank bottom. The same assumption was made (as in ref. 9) that all heat 
input below the point x = 13. 5 inches would be assigned to bottom heating. In a like 
manner, all heat input above the point x = 13, 5 inches was assigned to the side wall 
heating. 

The right hand side of equation (D5) then becomes 

13. 5 < Xg < L 
0 < Xg < 13. 5 

and similarily, the right hand side of equation (D6) becomes 


f Pj^dV = Qp(13.5) 

•'v 

= Qj^(Xs) 
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f dS = Qj^(Xg) - Qj^(13. 5) 13. 5 < < L 

•'S 

= 0 0 < < 13. 5 

^ s 

where Q^CXg) is obtained using equation (4). 

An average wall heat flux q^ required for the evaluation of the exponent n (see 
ref. 9) was then obtained from the relation 

Q^(L) - Qj^(13.5) 

"" s(L) - S(13. 5) 

where the surface area S was calculated from the tank geometry. 
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TABLE I. - FLOW TEST CONDITIONS 


^eat flux 
configuration 

Average wall 
heat flux, 
Btu/(sq ft) (sec) 

Average bottom 
heat flux, 
Btu/(sq ft) (sec) 

Initial height 
of liquid, 

L, 

in. 

Initial 

tempera- 

ture, 

'^in» 

°R 

Tank 

pres- 

sure, 

Prji, 

psia 

Satura- 

tion 

tempera- 

ture 

rise, 

ATg, 

°R 

Flow 

rate, 

w, 

Ib/sec 

Average 

specific 

heat, 

Btu/^b)(°R) 

Average 

density, 

P, 

Ib/cu ft 

Average 
heating 
rate ratio, 
Qj/Qh 

L - L 

2.81X10'® 

1. 10X10"^ 

44.5 

38.70 

33.39 

3.41 

0. 0405 

2.45 

4.32 

6. 840 

ML - L 

4.07 

.45 

45.1 

38.32 

34.40 

4.02 

.0417 

2.45 

4. 33 

2. 890 

MH - L 

6. 37 

.53 

46. 1 

38. 39 

33. 40 

3.95 

.0424 

2.48 

4.32 

1.519 

H - L 

7. 19 

.52 

44.7 

38.44 

34. 22 

3. 63 

.0411 

2. 50 

4.31 

1,204 

L - ML 

1.05 

3.61 

45.1 

38.28 

34.40 

4.06 

.0410 

2.45 

4. 33 

2. 617 

L - MH 

1. 22 

4, 94 

44.3 

38. 16 

34.31 

4. 16 

.0407 

2.45 

4.33 

1.613 

L - H 

1.74 

9.75 

44.7 

38.40 

34. 40 

3. 94 

.0415 

2. 50 

4. 30 

1. 577 

H - H 

8.23 

7.48 

44.2 

38.45 

34.40 

3.89 

.0416 

2. 53 

4. 29 

1.144 


^Average wall heat flux - average bottom heat flux; 


H, high; L, low; ML or MH, intermediate values. 



TABLE n. - TANK GEOMETRY 



= 4. 83 in. 
= 16 in. 

= 1. 41 in. 
X 2 = 9. 31 in. 
Xg = 20. 62 in. 


Region 

Equation of tank profile 

Geometry 

0 < X < x^ 

,w . [-x^ 

Sphere 

Xi < X < Xo 

r(x) = X + 2 

Cone 

1 — — £x 

Xg < X < Xg 

r(x) = j^-x^ + 2xgX + (r^ - x^J^^ 

Sphere 

Xg < X < L 

r(x) = r j 

Cylinder 


^Regions are distinguished by tangent points. 






TABLE m. 

- LOCATION OF TEMPERATURE PROBES 




Height from 


Probe location 




Height from 

Probe location 



bottom of tank. 







bottom of tank, 





' I Carbon resistor 


Thermocouple 



1 in. I Carbon resistor 

Thermocouple 


Axial 

Radial distance 

Axial 

Radial distance 

I Outer wall 

Radiant 

Axial 

Radial distance Axial 

Radial distance 

Outer wall 

Radiant 


from wall, in. 


from wall, in. 



heater 


from wall, in. 

from wall, in. 


heater 


0. 125 lo. 25 


0. 125 1 0. 25 





0. 125 1 0. 25 

0. 125 1 0. 25 






TABLE IV. 


NO- FLOW TEST CONDITIONS 


^eat flux 

Initial height 

Average 

Average heat 

Time period. 

configuration 

of liquid, 

heating rate 

input rate 

m’ 


L, 

of liquid, 

from heaters, 

sec 


in. 


Qh> 




W 

W 


L - L 

42.7 

35. 5 

4.0 

2290 

H - L 

42.3 

174.6 

132.0 

1155 

L - H 

41.7 

113.4 

68.2 

2000 

L - H 

26.0 

100.4 

64.4 

1740 

H - H 

42.0 

223.3 

195,0 

1680 

H - H 

24.9 

107.9 

100.2 

1680 


^Average wall heat flux - average bottom heat flux: H, high; L, low; 
ML or MH, intermediate values. 
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Figure 5. - Comparison of measured and predicted liquid hydrogen temperature histories at severai thermometer locations during outflow. 
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Heating error, (Qj - Q|-|)/Q| 



Figure 10. - Heat input rate distribution for 
boiloff test. 



Average heat input rate from heaters, Q|.j, W 


Figure 11. - Comparison of heating rate of liquid with heat input rate from heaters for 
no-flow tests. 


NASA- Langley, 1966 E-3139 



”The aeronautical and space activities of the United States shall he 
conducted so as to contribute ,,, to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof/* 

— ^National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washinston, D.C. 20546 



